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Abstract 

A detailed study of the effect of the ??-nucleus final state interac¬ 
tion (FSI) in the pd —> ^Her/ reaction close to threshold is presented. 

The FSI is incorporated through a T-matrix for r/^He elastic scatter¬ 
ing, constructed using few body equations. This T-matrix accounts 
for off-shell and binding effects in r^-nucleus scattering. The energy 
dependence of the data on the pd —> ^Hery reaction near threshold 
is reproduced only after including the FSI. The off-shell and binding 
effects in yy^He scattering are found to be important. Given the uncer¬ 
tainty in the knowledge of the elementary ry-nucleon interaction, the 
sensitivity of the pd —> ^Hery cross section to different prescriptions 
of the elementary t-matrix, trjN^riN is also discussed. 

PACs numbers: 13.75.-n, 25.40.Ve, 25.10.-l-s 
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1 Introduction 

Experimental data on the pd ^ ^Hep reaction close to threshold have re¬ 
vealed some surprising features [|^ Q . Inspite of the large momentum transfer 
involved in p production as compared to that in pion production, the cross 
section for pd ^ ^Hery is large and comparable with that for pd ^ ^Hevr^. 
The energy dependence of the two reactions near threshold is also very dif¬ 
ferent. The pd ^ ^Hep reaction shows a much rapid variation, with the 
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threshold amplitude falling by a factor of 3.75 over an rj centre of mass mo¬ 
mentum of 75 MeV/c. 

The observed features of the pd ^ ^He p reaction have been attributed 
to the strong hnal state ? 7 -nucleus interaction. This reaction was studied 
in Ref. where the hnal state interaction (FSI) was incorporated in an 
approximate way through an enhancement factor. The cross section was 
factorized in terms of the amplitude for the reaction p d —> ^He p with plane 
waves for the p^He in the hnal state and an S-wave FSI factor written in 
terms of the on-shell p-nucleon amplitude. The energy dependence of the 
cross section data was reproduced, though its predicted absolute value was 
lower by a factor of 2.5. In yet another work Q, the scattering length of the p 
meson on helium was calculated using multiple scattering theory, which was 
then used to calculate the FSI factor for the p d —>■ ^He p cross section. The 
FSI factor with an p-nucleon scattering length = (0.291, 0.36) leading to 
a^i3{-0.89, 1.8) was found to give a good ht to the pd ^Hep data. 

The pd ^ ^He p reaction at threshold and higher energies has also been 
studied theoretically to investigate the reaction mechanism involved ||, ||, 0]. 
Due to the large mass of the p meson, the momentum transfer involved in 
this reaction is large. It is 900 MeV/c at threshold and reduces to 500 MeV/c 
by about 1 GeV above threshold. As a result of this, it was shown in Ref. 
[^] that the three-body mechanism which allows the momentum transfer to 
be shared amongst three nucleons dominates. The one and two body mech¬ 
anisms were found to underestimate the experimental cross sections by more 
than two orders of magnitude. The three nucleons share the large momen¬ 
tum transfer through a two step process where the incident proton interacts 
with a nucleon in the deuteron to produce a pion which then interacts with 
the other nucleon in the deuteron to produce an p meson. The pd ^ ^He p 
reaction, thus, proceeds via the NN —>■ Tid and ttN pN reactions. Though 
the pd ^ ^Hep reaction can in principle proceed through the one and two 
body mechanism, considering the support in existing literature |^, Q for the 
three body mechanism and the fact that the low momentum components of 
the nuclear wave function are picked up in this way, at least at energies very 
close to threshold it seems justihed to consider this mechanism to be the only 
major reaction mechanism for the pd —>■ ^Hep reaction. 

In the present work we study the pd —> ^Hep reaction near threshold 
using the three body mechanism mentioned above. Our main objective is to 
investigate the FSI in this reaction in a rigorous way. The above mentioned 


2 


theoretical works in the literature either neglect the FSI or incorporate it in 
an approximate way using on-shell amplitudes. We express the rj ^He relative 
wave function in terms of the Lippmann Schwinger equation involving the 
T-matrix for elastic scattering. This T-matrix is evaluated using a 

method of few body equations 1 , 0,0 which will be described in detail in 
the next section. 

In section 2 we present the formalism for the calculation of the pd —> 
^He Tj cross section including the FSI. The production mechanism for pd ^ 
^Hep reaction within a two-step model is described in section 3. In section 
4 we present the results. We reproduce the shape and magnitude of the 
experimentally measured scattering amplitude. The off-shell effects in p^He 
scattering are found to be important in producing the energy dependence of 
the cross section. 


2 Final state interaction 


The transition matrix for the reaction pd ^ ^Hep, which includes the in¬ 
teraction between the r] meson and ^He is given by, 

T = < ih) ; I Tpd^SHerj | fcp ; mi 7712 > (1) 

where mi, m 2 and m 3 are the spin projections of the proton, deuteron and 
helium respectively, kp and kp are the momenta of the particles in the initial 
and hnal states. The hnal state p ^He wave function consists of a plane 

wave and a scattered wave, and can be written as. 


< \ -< kp \ + 


dq < kp I Tp3He I q > 
(27r)3 E{kp) - E{q) + ie 


( 2 ) 


where Tp3He is the T-matrix for p ^He elastic scattering. Replacing the above 
wave function in Eq. (P, we get. 


T 


kp , m3 I Tpf^_^ 3 j^^p I kp , mi m2 ^ T 


( 3 ) 



dq < kp; ms\ Tp3He | 9 ; ^3 > 
(27r)3 E{kp) - E{q) + ie 


< q , mg \Tpd^^Heri 



mi m2 > 
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The matrix elements < \Tpd^iHe'n\ > in the above equation correspond to 
the Born amplitude for the pd ^ reaction. We calculate these matrix 

elements using a two step model which will be discussed in the next section. 

The T-matrix, in Eq. (|^) for elastic scattering is evaluated 

using four particle equations for the p{3N) system. For practical conve¬ 
nience, the evaluation of is done within a Finite Rank Approximation 

(FRA) approach. This means that the nucleus in the elastic meson-nucleus 
scattering is always in its ground state. The shortcomings of the FRA for 
the 77 -deuteron system have been investigated in Ref. [|T^]. However, for the 
low energies concerned in this work and also as mentioned in Ref. [|^, it 
seems justihed to use the FRA for the rj — ^He system. We write the target 
Hamiltonian Ha as, 

Ha ~ e\^po >< V’ol (4) 

where V’o is the nuclear ground state wave function and £ the binding energy. 

Within this approximation, the p ^He T-matrix is given as 8001. 


T{k', k ; z) 


< fc'; V’o I T^{z) \ k]'il)o > -F 

f dk" < k'-,tl)o\ T%z) I Id '; i/jo > 


(2vr) 




fc //2 . 


T{k", k ; z) 


2fi 


2/i 


(5) 


where z = E — \e\+ E is the energy associated with //-nucleus relative 
motion and fi is the reduced mass of the //-nucleus system. The operator 
T° describes the scattering of // meson from nucleons hxed in their space 
position within the nucleus. The matrix elements for T° are given as. 


< fc'; -00 I T^{z) \ k;^/Jo >= j dr\ T° (fc', k]r]z) ( 6 ) 


where. 


T\k',k-f-z)=Y.T^{k',k-T--z) 


(7) 


2 = 1 

is the t-matrix for the scattering of the //-meson from the nucleon in 
the nucleus, with the rescattering from the other (A- 1 ) nucleons included. It 
is given as. 


Tf {k', k]fi]z) = ti{k', k]fi]z) + 


dk" ti{k', k" ]fi] z) 


( 8 ) 
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The t-matrix for elementary r^-nucleon scattering, ti, is written in terms of 
the two body rjN matrix as, 

U{k', k]fi]z) = tr,N{kk k]z) exp[i{ k-k') ■ n] (9) 


The ^He wave function - 00 , required in the calculation of is taken 

to be of the Gaussian form. 

Since there exists a lot of uncertainty in the knowledge of the 77 -nucleon 
interaction, we use three different prescriptions of the 77-N t-matrix, 
leading to different values of the 77-N scattering length. We give a brief de¬ 
scription of these models oftrjN^rjN below. The elementary matrix, n^tjN 
is then used in the evaluation of the T-matrix for rj ^He scattering. In Ref. 

a coupled channel t-matrix including the vrN and 77N channels with the 
Sii - 77N interaction playing a dominant role was constructed. The t-matrix 
thus consisted of the meson - N* vertices and the N* propagator as given 
below: 


t 


riN^riN' 


k', k] z) 


(fc/2+^2)^iV*Gl (P + /32) 


( 10 ) 


with, 

V* (^) = ( 2 ; - Mo - E^(z) - Erj(z) + ie (11) 


where Sa(^) (a = vr, 77) are the self energy contributions from ttN and rjN 
loops. The parameters of this model were chosen to produce an 77 N scattering 
length, a^jjv = (0.75, 0.27) fm which is in agreement with some other analyses 
available in literature |1l4|, |T^. We will refer to this parameter set as ‘Fix et 
al. (I)’ while discussing the results. 

In Ref. the authors also perform calculations with the choice of 

meson-nucleon cut-off parameters as in Ref. [M. These parameters lead to 


a»? 7 v = (0.88, 0.41) fm. We refer to this parameter set as ‘Fix et al. (II)’. We 
present our calculations for the pd —>■ ^B.ep reaction with the t-matrix of 
Ref. 1^ using the two different sets of paramteres mentioned above. 

We also present results using the 77 -N t-matrix of Ref. |]^ which gives a 


much smaller value of a^^Ar, namely, a^Af = (0.28, 0.19) fm. In this model the 
vrN, 77 N and ttA channels were treated in a coupled channels formalism. It is 
important to note that the parameters in all the coupled channel t-matrices 
mentioned above are adjusted to reproduce the data on the N p N 
reaction, but the values of a^jN deduced from them are different. 
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3 Production mechanism 


As mentioned in the Introduction, we assume the r] production in p d —^He rj 
to proceed through a two-step process via the NN —> vrd and ttN —>■ t]N 
reactions as shown in Fig. 1. The amplitude for the pd ^ ^Hep reaction 
appearing in Eq. (^) can be written within this model as, 

< \Tpd-^3Herj\ >= i [ <pn\d><TT+d\Tpp^^+d\pp>il2) 

J [ZTT) [ZTT) 

1 , , . s , 

, . . <pp\ T^N^rip tt ^ > < He\pd> 
-ml + %e) 

where the sum runs over the spin projections of the intermediate particles. 



Figure 1: Diagram of r] production in the pd ^ ^Her^ reaction with a two 
step process. The ellipse indicates the hnal state interaction of ^He and r]. 

The spin projections and momenta of the interacting particles are as shown in 
Fig. 1. /Cjr is the four momentum of the intermediate pion which could either 
be or 7r°. In the case of an intermediate 7r°, the matrix element for pd —> 
^He r] is half of that written above for 7r+. Hence we calculate the T-matrix as 
in Eq. ( [T^) and multiply it by a factor of 3/2 to account for the intermediate 
7r°. Each of the individual matrix elements in the above equation is expressed 
in terms of partial wave expansions. The matrix elements for the pp —>■ 7r+d 
reaction, parametrized in terms of the available experimental data are taken 
from Ref. [0. For the vr+n —pp reaction, we use the coupled channel t- 
matrix of Ref. 10. mentioned in the previous section. The matrix elements 
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< pn\d > and <^}le\pd > consist of the denteron and helinni wave fnnctions 
in momentnm space. We nse the denteron wave fnnction from Ref. [W 
where an analytical parametrization of it was done with a Paris potential. 
This wave fnnction reprodnces the known low energy properties and the 
electromagnetic form factor of the denteron well. For the ^He wave fnnction, 
we nse the parametrization given in Ref. pO • The valnes of the parameters in 


were obtained by htting the wave fnnction to the variational calcnlations 
in using the Urbana force. The details of Eq. (|^ are 


of Schiavilla et al. 


given in the appendix. 


4 Results and Discussion 





Figure 2: The square of the pd ^ ^Her/ amplitude dehned in Eq. as a 
function of the r] momentum in the centre of mass. The data is from Refs 
[0, |[. The dotted line is the calculation of the present work without in¬ 
cluding the FSI. The solid, dash-dotted and dashed lines are the calculations 
including FSI with different prescriptions of elementary t-matrices. 

The reaction pd —> ^Hep has been studied at Saturne [^, ^ for proton 
energies between 0.2 and 11 MeV above threshold. Taking out the phase 
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space factor, the spin averaged amplitude can be defined as, 

I ^ |2 _ fcp da 

dkl cm 


( 13 ) 


where kp and kp are the proton and rj momenta in the centre of mass system. 
The data on |/p (see Fig. 2) drops rapidly (by about a factor of 3.75 ) from 
threshold to 0.4 fm~^ momentum (corresponding to 11 MeV energy) above 
threshold. In Fig. 2 we compare our calculations of |/p (at 6p = 180°) with 
and without the inclusion of the hnal state interaction (FSI), with the 
data from Refs ii- The dotted line in Fig. 2 is our calculation without FSI 
and can be seen to be a constant as a function of energy. The dash-dotted, 
solid and dashed curves are the results obtained using the t-matrix of Ref. 
r^ , with parameter sets (I) and (II), and that of Ref. 0 respectively. We 


see that the FSI is responsible for changing the shape of |/p from a constant 
to a rapidly falling one as a function of energy. 

In what follows, we shall study the off shell and the binding energy effects 
in the FSI. We also calculate the angular distribution and total cross-sections 
for the p d —> °He rj reaction. In all these calculations we use the prescription 
of Ref. [0 with parameter set (II) for the elementary t-matrix, tpN^rjN- 
We arrive at the same results qualitatively, in all the calculations if we choose 
the parameter set (I) of the same t-matrix or the elementary matrix, tpN^rjN 
of Ref. 0. 

In Fig. 3 we study the off-shell effects in the FSI. As seen in section 2, 
we describe the p^He hnal state interaction through an off-shell T-matrix 
for r] ^He scattering. Previous estimates of FSI in literature have been made 
using on-shell amplitudes and hence it is important to check the validity of 
such an approximation. The p°He T-matrix appears in an integral in Eq. 
(^ which can be split into the principal value and pole term. Retaining 
only the pole term in Eq. (|^) and setting the principal value (which involves 
the off-shell p°He scattering) to zero, we get the dash dotted line in Fig. 3. 
The pole term alone is unable to reproduce the shape of the data and in 
fact reduces the magnitude of the results obtained without FSI. We hnd that 
tpN^riN of Ref. 10, with parameter set (I) and (II) and that of Ref. 0 


have similar on-shell behaviour meaning that the same dash dotted curve of 
Fig. 3 represents the pole term calculation for the 3 different prescriptions of 
The solid line represents the full calculation using the prescription 
of Ref. 01 with parameter set (II) for the elementary t-matrix. Thus we see 
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p d ^ He p 


Mayer et al. 
Berger et al. 
Tpsi(full) 



0.1 


0.2 

Pn (fm ') 


0.3 


0.4 


Figure 3: The solid and dotted line and the data are as in Fig. 2. The dash 
dotted line is the FSI calculation retaining only the pole term in Eq. (j^). 
We have used the elementary t-matrix of Ref. with parameter set (II). 


that including the off shell effects produces the proper energy dependence of 
I/P in addition to increasing its magnitude as compared to the pole term 
calculation. The other two prescriptions of (Fix et al. (I) and 

Bhalerao et al.) show similar off-shell effects. However, as can be seen from 
the curves in Fig. 2, the increase in magnitude as compared to the pole term 
is smaller. 

Next, in Fig. 4, we study the effect of the FSI on the angular distributions 
at different energies. As observed in Fig. 2 too, we see that the FSI increases 
the magnitude of the cross sections with the increase being maximum at 
threshold (Fig. 4a). The angular dependence of the cross sections without 
FSI (dashed lines) is isotropic at threshold and deviates by a small amount 
from the isotropy with increase in beam energy. The small anisotropy at a few 
MeV above threshold is somewhat amplihed by the hnal state interaction. 
Thus, as one goes to higher energies (Fig. 4c), though the FSI does not 
increase the magnitude of the cross section too much, it does change the 
angular dependence. 

In Fig. 5 we show the angle integrated total cross section as a function of 
beam energy. The total cross sections without FSI are shown by the dashed 
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Figure 4: Angular distribution of rj in the pd ^ ^He?) reaction at different 
beam energies. The solid (dashed) curves are the calculations with (without) 
FSI included. 

line. The solid line is obtained by numerically integrating the angular dis¬ 
tribution with FSI included. The dash dotted curve indicates the total cross 
section calculated with FSI included but assuming the angular distribution 
to be isotropic (i.e. atot = 4:7i{da/dfl)g^=igQo). The data of Ref. [Q] seems 
to indicate a negligible forward-backward asymmetry (which is consistent 
with zero within 5% at all energies from 0 to 11 MeV above threshold) in 
the pd —> ^Hep reaction. Within the two-step FSI model of the present 
work, however, we hnd small deviations from isotropy at energies away from 
threshold. Since the total cross section assuming isotropic angular distribu¬ 
tion gives better agreement with data, it seems that the two step model of 
the present work underestimates the cross section at forward angles. The 
data on the squared amplitude | / p at 180° is however very well produced 
as seen earlier in Fig. 2. 

Finally, in Fig. 6 we study the effects on the pd —> ^Her/ reaction, of 


10 








0.8 


0.6 




0.4 


0.2 


p d ^ ‘ He T] 


• Mayer et al. 

-NoFSI 

-a = 4jr(do/d£l)jgg“ 

- a = J(da/d£2) 


--f- 



892 


894 896 898 

Tp (MeV) 


900 


902 


Figure 5: Total cross section for pci —> ^Hep as a function of beam energy. 
Solid curve is obtained by numerically integrating the calculated da/dVL with 
FSI over all angles. The dash dotted curve is obtained by multiplying da/dfl 
(with FSI) at 180° with dvr, thus assuming the angular distribution to be 
isotropic. The dashed curve is the calculation without FSI. The data is from 
Ref. i. 


including the FSI which incorporates the binding effects in °He. Once again 
we plot the amplitude squared |/p calculated for three different elementary 
t-matrices, The solid curves are the calculations using e = -7.718 

MeV in Eq. (|^) for T^SHe- The dashed curves represent the calculation with 
£ = 0, i.e. no binding effects in the hnal state interaction. The calculated 
results without the binding effects are larger than those which include them. 
The magnitude of the increase, however, depends upon the elementary t- 
matrix, tjjN^riN, under consideration. It is nearly a factor of 2 when we use 
the t-matrix of Ref. 1^, with parameter set (I), a factor of 1.25 when we 


use the same t-matrix with parameter set (II) and negligibly small for the 
t-matrix of Ref. ITTI . 
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Figure 6: Effects of including the binding energy of ^He in the calculation of 
FSI between ^He and rj. The solid and dashed curves are calculations includ¬ 
ing FSI with the binding energy e = -7.718 and zero respectively (see Eq.(§)). 
The effect is shown for three different elementary t-matrices, trjN^rjN- 

5 Summary 

The pd ^ ^He p reaction has been studied in the present work within a two 
step model, incorporating the hnal state interaction (FSI) of the ^He nucleus 
and p meson in a rigorous way. The peculiar behaviour of the cross section 
for this p producing reaction as compared to a similar pion production re¬ 
action pd ^ ^He7r° is seen to originate due to the interaction between ^He 
and p. The FSI changes the energy dependence of the squared amplitude 
from a constant (without FSI) to one which falls by a factor of about 4, from 
threshold to 11 MeV above threshold. We incorporate the FSI through an 
off-shell T-matrix for p ^He elastic scattering. This T-matrix is evaluated by 
numerically solving few body equations which include the nuclear binding 


12 







effects. Both the off-shell as well binding effects in rj ^He scattering are found 
to be important in the calculation of the pd ^ ^He p reaction near thresh¬ 
old. Earlier investigations of this reaction involving on-shell and approximate 
ways of calculating the ESI should hence be treated with caution. 

The p nucleus interaction has generated a lot of interest in the past few 
years, particularly due to the possibility of forming p mesic nuclei. Since it 
is difficult to obtain data on elementary p nucleon scattering, little is known 
about the pN interaction. The scattering length in pN scattering is a much 
debated quantity and different estimates and limiting values (for the possible 
formation of an p mesic nucleus) of this parameter exist in literature. Within 
the models used in the present work we hnd that values of the pN scattering 
length. Re ar,N ~ 0.75 to 0.9 and Im ar,N ~ 0.3 to 0.4 lead to a good 
reproduction of the pd ^ ^He p data near threshold. 
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Appendix 

In what follows, we discuss in detail the constituents of the Born ampli¬ 
tude (Eq. (|l^)) for the pd —>■ ^Hep reaction. To start with, we write the 
deuteron wave function in Eq. ([T^) as, 

< pn I d > = ^ I XI < 1/2 mp 1/2 11 m2 > (A.l) 

V ^ mn A/dVr 

+ X < 1/2 mpl/2 mn\lMs><l Ms 2 mi\l m 2 > p 2 ,mi(Pi) (f> 2 iPi) } 

mi 

where m^ and rup are the spin projections of the off-shell neutron and proton 
respectively in the intermediate state and m 2 is the spin projection of the 
target deuteron. The factor A comes from isospin overlap. 0f(pi) is the 
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deuteron wave function in the partial wave and pi is the relative momen¬ 
tum of the p-n pair inside the deuteron. Writing 


00 (Pi) 


n 

(2/x)i E 


i=i 


Pi + mj 


(A. 2 ) 


02 (Pl) 


n 


(2/ir)’ E 


9 I 2 ’ 

Pi + rrij 


(A.3) 


the parameters Cj , Dj and rrij for the Paris potential are given in Ref. 
The ^He wave function is written as, 


IP 


< ^He \ pd >= < 1 m 2 1/2 m^ 1 1/2 m 3 > -|- (A.4) 

< Im'a 1/2 m'p \ 3/2 m >< 2 m^ 3/2 m 11/2 m 3 > X 2 (P 2 ) l 2 ,m;(p 0 ) 


where m 2 and m^ are the spin projections of the off-shell deuteron and proton 
respectively. The spin projection of ^He is m 3 as shown in Fig. 1. Xi{P 2 ) 
is the helium wave function in the partial wave and p 2 is the relative 
momentum of the p-d pair inside ^He. 


Xl{P2) 


E 


P2 + mf 


(A.5) 


The parameters a* and mj given in Ref. are chosen corresponding to p-d 
clustering in ^He. The normalization of the wave function is such that, 


J plip2 {xofe)'^ + ,\ 2 (P 2 )^} = 1.5 


(A. 6 ) 


The four momentum k,r = { E^r , /Ctt } appearing in the pion propagator 
in Eq. (^) is written using energy and momentum conservation at the 
n ^ pp and pp —> vr+d vertices respectively. Thus, 

Ett = + '^Ehc ~ -^Efi (A.7) 


and 



2r - - 

-kr,+pi+P2 


(A. 8 ) 
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